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ABSTRACT 



Context. Despite the privileged position that Lyman-o' (Lya) emission line holds in the exploration of the distant universe and modern 
observational cosmology, the origin of the observed diversity of Lya profiles remains to be thoroughly explained. Observations of 
nearby star forming galaxies bring their batch of apparent contradictions between Lya emission and their physical parameters, and 
call for a detailed understanding of the physical processes at work. IZw 18, one of the most metal-poor galaxies known is of particular 
\^ . interest in this context. 

Aims. Fitting the Lya spectrum of IZw 18 to understand the origin of the damped absorption profile and its spatial variations across 
the NW region. To establish a general picture of the physical parameters governing the Lya strength and profile both in local and in 
high-z galaxies. 

^ , Methods. We use a 3D Lya radiation transfer code to model Hubble Space Telescope (HST) observations of IZw 18. Observational 

' constraints of relevant parameters such as dust or Hi column density are derived from previous studies and from the present analysis. 

' Different geometrical configurations of the source and the neutral gas are explored. 

\ Results. The integrated Lya profile of NW region of IZw 1 8 is reproduced using the observed small amount of dust (E{B - V) « 0.05) 

and a spherical Hi shell with A'hi = 6.5 x 10^' cm"-. Such a high column density makes it possible to transform a strong Lya emission 
(EWlju = 60 A) into a damped absorption even with a small extinction. When a slab geometry is applied and a given line of sight 
is chosen, the Lya profile can be successfully reproduced with no dust at all and A'hi = 3 x 10^' cm"-. The spatial variations of 
, the profile shape are naturally explained by radiation transfer effects, i.e. by scattering of Lya photons, when the observed surface 

brightness profile of the source is taken into account. In the case of outflowing Inter Stellar Medium (ISM), as commonly observed in 
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Lyman Break Galaxies (LBGs), a high A'hi and dust content are required to observe Lya in absorption. For nearly static neutral gas as 
observed in IZw 1 8 and other low luminosity galaxies only a small amount of dust is required provided a sufficiently high A'hi covers 
the galaxy. We also show how geometrical and apertures effects affect the Lya profile. 

Key words. Galaxies: starburst - Galaxies: ISM - Ultraviolet: galaxies ~ Radiative transfer - Galaxies: individual: IZw 18 

1 . Introduction low redshift (z ~ 0.2-0.35 IPeharveng et ani2008h . thanks to its 

wide FOV and spectroscopic capabilities in the UV. 

The detection of high redshift (z) galaxies has become, through Earlier studies of nearby galaxies using mostly the UV ca- 

^ the last decade, a routine fact, although the discovery of primeval pabilities of lUE and HST recognized quickly that Lya emis- 

galaxies that are forming their first stars still remains a challenge. ^^^^ ^^s fainter than naively expected from recombination the- 

Depending on the selection techniques, mainly two classes of ^^^j t^at the Lya hne showed a great diversity of pro- 

galaxies are found: Lyman Break Galaxies (LBGs) selected by fijgg f^om absorption to emission (e.g. Meier & Terlevich" 198 Ij 

their UV continuum break, and Lyman-Alpha Emitters (LAEs) Hartmann et al. 1984; Deharvens et al. 1985; Hartmann et al 

selected upon their Lya emission line. The situation was, how- I988c ITerlevich et al.1 11993; Le aueux et al.l 119951: IThuan et al. 

everfar different before this successful era. Partridge & Peebles ^ggj^ iThuan & Izotoviri997: .Kunth et al.lll998ir Later. HST has 
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•112671) were the first to estimate that young distant galax- allowed one to map the spatial distribution of Lya emission and 

les should be detectable through a strong Lya emission, absorption, of the stellar sources, nebular emission, and dust 

Nevertheless it took nearly thirty years until such populations (^gg Kunth et al. 2003; Mas-Hesse et al. 2003; Hayes et al. 20051 

could be found at z ~ 2-7, thanks in particular to instruments 2008 HAtek et alEoOSh . Despite these information, no clear pic- 

with large fieldof view (FOV) and 4-8m class telescopes (see j^j^g ^as yet emerged explaining consistently the Lya and related 



e.g. |Huetalj[T998; Kudritzki et al. 2000; Malhotra & Rhoads observations in nearby starbursts 

2002liAiikietal. 2003; TaniHUchietal. 2005; Shimasakuet ah ^^^^^^ ^^jl ^^^^^ theoretically and 

2006; Kashikawaetal. 2006; Tapkenetal. 2006; GronwalletaU g ideally that different physical processes affect the Lya 

2007t lOuchietalJ 12008; NilssonetalJ |2008). Only recently, -J^^- ^^^^jg and "visibility" (i.e. detection fre- 



GALEX has provided for the first time a comparable survey at ^^^g^^y ^^^^^ starburstsfl: destruction of Ua photons 
by dust (cf. Neufe ldj 1 19901; fchaii ot & Fall Il993h . veloc- 



Based on observations made with the Hubble Space Telescope ob- 



tained from the ESO/ST-ECF Science Archive Facility ' See lSchaererl ( 120071) for an overview. 
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ity fi elds in the ISM (iKunth et all 119981: ITequeux et al 
1995"), an inhomogeneous ISM ('Neufeld"199l1 :lGiavalisco et al 
1996; Hansen & Peng Oh 2006), underlying stellar absorp- 
tion (IVa lls-Gabaud 1993) , and star formati on duty cycles 
(IValls^abaudl993^,Malhotra & Rhoadsll2002h . 

Also an "unifying" scenario to explain the observed diver- 
sity of Lya profiles in terms of an evolutionary sequence of 
starburst driven super-she lls/superwind has been presented by 
iTenorio-Tagle et aTl (Il999l)_ and has been confronted with local 



Table 1. HST observ ations of IZw 18. R eferen ces: (1) 
iMas-Hesse et al.' p003h : (2) iBrown et al.1 ([2002h : (3) 
Cannon et al. (2002|) 



starburst observations b y 'Mas-Hesse et al.' ("2 003). For distan t 
galaxies. ISchaerer & Ve rhamme (2008): Verhar nrne et al.l (|2008|) 
have recently shown - using radiation transfer models and em- 
pirical constraints - that Lya line profiles of high-z LBGs and 
LAEs can well be reproduced and that the diversity of Lya 
from emission to absorption is mainly due to an increase of the 
dust content and the Hi column density. Despite this progress, a 
global picture identifying the main processes and explaining this 
diversity also in a quantitative manner is still missing for Lya in 
local/nearby galaxies. Furthermore differences between the high 
and low redshift samples - if any - must be understood. We here 
provide a first step towards these goals by examining and model- 
ing one of the most metal-poor star forming galaxies in the local 
Universe, IZw 18, and by putting i t into context. 

Since its discovery bv lZwickvl(ll966l) . IZw 18 has been stud- 
ied extensively, and it remains one of the most metal-poor galax- 
ies kn own today (Skillman & Kennicutt 1993; Izotov & ThuanI 
119991) . Its main Hii region (called the NW region, cf. Fig. 
[Ij showing strong optical emission lines is clearly a site of 
very recent (< 10 Myr) and/or ongo i ng massive st a r for- 
mation (cf. iHunter & ThronsonI 119951: Ide Mello et all Il998t 
IBrown et all I2002h . Ther efore the find i ng of abroad damped 
Lyor absorption fine by iKunth et all (11994 Il998l) came as 
a surprise, where strong emission was predicted, given the 
strong optical H recombination lines and t he low dust con- 
tent expected for suc h low metallicities (cf. iKunth et al.l [79941: 
lTerlevichet"aLlll993l) . Observations of SBS 0335-052, nearly as 
metal-poor as IZw 18 , showed later a similarly broad profile 
dThuan & Izotoviri997l) . However, since compared to IZw 18, 
SBS 0335-052 has a higher extinction and is now known to har- 
bour more dust b oth in absolute terms (dust mass) and i n relative 
terms (L,R /L.rv) dThuan et alJl999l: lHo" uck et al.l2004t lWu et alj 



120071: lEngelbracht et al.ll2008l) . it is a priori more challenging to 
explain Lya absorption in IZw 18 than in SBS 0335-052. For 
these reasons IZw 18 represents an ideal test case to understand 
how strong intrinsic Lya emission is transformed to the observed 
broad Lya absorption, in a dust poor (but not dust-free), very 
me tal-poor galaxy. 

iKunth et ai] j l994 11998 ') suggested that both dust absorp- 
tion and multiple scattering of Lya photons out of their narrow 
(2.0"x2.0") GHRS/HST aperture as the most natural explana- 
tion for the observed weakness of Lya in IZw 18. They also 
noted that all galaxies showing Lya absorption (4/8 in their small 
sample) showed nearly static neutral gas, which must increase 
the mean free path of Lya photons. However, the lUE spectrum 
of IZw 18 taken with an entrance hole of 20"xl0" shows ba- 
sically the same profile, indicating that at least over 5-10 times 
larger scales no significant amount of Lya emission is recov- 
ered. In any case, whether quantitatively these explanations are 
viable remains to be seen. This is one of the concrete goals of 
the present paper. 

To address the above questions we will use the most re- 
cent observations of IZw 18 and our state-of-the-art 3D Lya and 
UV continuum radiation transfer code MCLya ( Ver hamme et al.l 
I2OO61) . This will in particular allow us to carefully examine in 



Instrument 


Filter/ 


Band 


Exposure 


Proposal 


Ref. 




Grating 




time [s] 


ID 




STIS 


G140M 


Lyo- 


1764 


GO-8302 


1 


STIS 


G140L 


Lyo- 


40360^ 


GO-9054 


2 


STIS 


F25SRF2 


FUV 


533 1~ 


GO-9054 


2 


STIS 


F25QTZ 


NUV 


5786 


GO-9054 


2 


WFPC2 


F487N 


H/J 


2500 


GO-6536 


3 


WFPC2 


F658N 


Ha 


4600 


GO-5434 


3 


WFPC2 


F450W 


B 


4600 


GO-5434 


3 


WFPC2 


F675W 


R 


2000 


GO-5434 


3 



" Total integration time in the 7 slit positions 



a quantitative manner the possible explanations leading to Lya 
absorption in IZw 18. Finally we will also discuss other nearby 
starbursts with Lya absorption, and place the local objects in a 
broader context. 

Our paper is structured as follows. In Sect. |2] we describe the 
main observations from HST and other facilities and summa- 
rize the main observational constraints. In Sect. [3] we set out to 
explain the Lya absorption in IZw 18, discussing our radiation 
transfer modeling tool, geometrical effects, and presenting mod- 
eling results for different ISM geometries. Our results for IZw 1 8 
are discussed and compared to other nearby and high-z starbursts 
in Sect.m In Sect.|5]we summarize our main conclusions. 



2. Observations 

The main observational data used in this paper are summarized 
in Tabled 



2.1. Spectroscopy 

We use in this work archival spectroscopic observations ob- 
tained, with the Space Telescope Imaging Spectrograph (STIS) 
onboard HST, under program GO-9054, by Brown et al. (2002). 
G140L grating was used combined with the 52"x 0.5" slit. IZw 
18 was spatially covered with seven adjacent slit positions along 
its main axis (see Fig.[Tli. Standard calibrations were performed 
using the CALSTIS pipeline (Ver 2.26), and exposures (two) for 
each position are registered and co-added. In addition, data were 
coiTected for geocoronal Lya emission by fitting and subtracting 
the nearby background regions in individual spectra. This cali- 
bration and spectra extraction were performed using IRAF and 
IDL routines. 

In Fig. 12 we show the spatial variations of the Lya profile 
across the NW region. Spectra were extracted from the seven 
adjacent positions of the STIS long slit covering the galaxy in 
the NE-SW axis, providing spatial information in two directions. 
Flux was then integrated in each slit along a 4 " aperture cen- 
tered on the NW component of IZw 18. Finally, an integrated 
spectrum of the NW region was also constructed from these in- 
tegrated slit spectra. The strength of the Lya absorption in these 
spectra is quantified by its equivalent width and corresponding 
Hi column density, A^hi, determined assuming a Voigt profile and 
b -2Q km s"' (cf. below). These values are reported in Table|2] 
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Fig. 1. FUV Image of IZw 18 showing the different observation 
apertures. The FOV is 16"x 18". The first slit position is in the 
NE direction and the seventh toward the SW. The integration is 
performed within a region of 4 " along the spatial axis of the 
North- West region and is marked with dashed lines. 
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Fig. 3. IZw 18 spectroscopic data. The figure presents a com- 
pilation of spectroscopic informations available for IZw 18. 
Geocoronal Lya emission has been subtracted from all the spec- 
tra. Dark solid line represents the best STIS spectrum around 
the Lya absorption. The red dotted line is the GHRS spectrum 
which covers a part of the absorption and a part of the UV con- 
tinuum. The blue dot-dashed line shows the STIS spectrum with 
a large wavelength coverage ([1 100-1750 A]) extracted from the 
center of the NW region. The green long dot-dashed line is the 
result of an integration over all the slits in the NW region, pre- 
viously shown in Fig. |2] It has been used to fit the UV contin- 
uum (black dashed line). All the spectra were then normalized 
to match the continuum value around 1280 A (i.e. at v ~ 16000 
km s '). References and legend are given in the inset. 
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Fig. 2. STIS spectra of the NW region of IZw 18 taken at dif- 
ferent locations (cf. Fig. [TJ. The slit positions 4 to 7 show the 
variation of the profile shape from the center to the edge of NW 
region. For comparison, the integrated spectrum in the seven slits 
over the center of NW region is over-plotted. All spectra were 
normalized to the continuum value determined in Fig.|3] For the 
sake of clarity, a cut-off is applied on geocoronal emission resid- 
uals and spectra are smoothed with a 3-pixels boxcar. 

Earlier, Lya o bservations of I Zw 1 8 were obtained by 
iKunth et al] (11994 and later on by iKunth et al] (1998) using 
the Large Science Aperture (LSA, 2"x 2") of GHRS onboard 



HST (see Fig. [B- 'Mas-HesseetaD (l2003h (hereafter MH03) re- 
observed the galaxy with better settings using STIS with G140M 
grating through a 52"x 0.5"longslit, translating to a spectral res- 
olution around 0. 15A (37 km s"' at Lya wavelength). The longer 
wavelength range of the STIS observations allows a better cov- 
erage of the Lya absorption red wing, as compared to GHRS 
spectrum, and confirms the large damped Lye absorption. In Fig. 
|3] we plot together the STIS and GHRS spectra. Due to differ- 
ences in the instrument apertures, the spectra had to be matched. 
To fit and estimate the UV continuum (dashed line) we used the 
archival STIS observations that allows a broad wavelength cov- 
erage to include the absorption wings. All the spectra were then 
normalized to the value where the Lya red wing reaches this 
continuum (~ 1300A). The different Lya profiles obtained are in 
good agreement. A correct estimation of the continuum around 
Lya is particularly important for the modeling of the Lya spec- 
tral profile (see Sect. 13.2b 

Clearly, Lya shows a broad absorption over the entire extent 
of the NW region. The width of the profile corresponds to an 
Hi column density of A^hi~ (0.3 - 3) x 10^' cm"^, in agreement 
with earlier determin ations (A^hi~ (1-0 - 3.2) x 10^' cm~^ from 
UV observations by 'Kunth et"an (1199 8'). although this method 
does not give systematically the true value of A^hi, as we will 
see later on. Beyond the scale of the Hii NW region (~ 250 pc), 
the UV-optical part of IZw 18 is known to be e mbedded in a 
large neutral Hi cloud extending over several kpc dvan Zee et al.l 
11998). Furthermore, the strength of the Lya absorption decreases 
clearly from the center to the border of the NW, as shown in 
Fig. 12] for slits 4 to 7 (slits 3 tol show also a slight decrease), 
corresponding to an apparent change of A^hi by up to a factor of 
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Table 2. Spatial variations of Lya properties in STIS slit posi- 
tions. Columns (2) and (4) indicate respectively the Hi column 
density and the Lya line equivalent width derived by fitting each 
Lya absorption with a Voigt profile with b - 2Q km s"'. The 
errors in cols. (3) and (5) are determined from the lower and up- 
per limits of the fits. Last column is the integrated flux in FUV 
(1500 A) image over the NW region and in simulated slits in 
order to match the aperture used for the extraction of STIS spec- 
tra. Same quantities are also given for integrated spectrum in the 
entire NW region and for the MH03 STIS spectrum. 
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~ 10. This systematic change of the Lya will be explained below 
as due to simple radiation transfer effects (Sect. [JX4b . 

Because of the configuration of their apertures, centered on 
the bright UV peak of the NW region, the GHRS and STIS 
MH03 spectra are in good agreement (cf. Fig. 13). This is also true 
for the central slit of the 2002 STIS data. The integrated spec- 
trum of the NW region shows residual emission at the center that 
may be due to the contribution of the external slits that show such 
emission and to geocoronal Lya residuals, as the STIS MH03 
and GHRS spectra have a better resolution which allow a more 
reliable correction. Since the profile shape in the wings remains 
the same, we will take the STIS MH03 spectrum as a proxy for 
an integrated spectrum of the NW region in our modeling. We 
derived for this spectrum an Hi column density of ~ 2 x 10^' 
cm-^ from Voigt profile fitting (cf. Table. |2|. 



2.2. Imaging 

UV images, part of the same observing program GO-9054 as that 
of the STIS spectroscopy, were retrieved from the ESO/ST-ECF 
archive. IZw 18 was observed with F25SRF2 filter with band- 
pass centered at 1457 A. Standard calibrations were performed 
through CALSTIS pipeline. Images were then corrected for mis- 
alignment, divided by the exposure time, and co-added. The final 
FUV image was then multiplied by PHOTFLAM and PHOTB W 
header keywords to obtain flux calibrated image. 

We also retrieved from the archive HST optical images ob- 
tained with the Wide Field Planetary Camera 2 (WFPC2) un- 
der programs GO-6536 and GO-5434. Data consist of Ha and 
Hj6 narrow band imaging and corresponding broad band contin- 
uum images (see Table[T]l. Data were first processed through the 
standard HST pipeline that gives images in units of counts per 
second. Multiplying by PHOTFLAM keyword gives fluxes in 
erg S-' cm-^ A"'. Finally all images are rotated and aligned to 
the same orientation and co-added in each filter using inverse 
variance weighting. OWe estimated the line flux contribution 



to the continuum images using the appropriate filter through- 
put ratios at Ho' and Hy6 wavelengths and filter width given by 
PHOTB W. Continuum images were scaled and subtracted from 
online images, then multiplied by the filter bandwidth to ob- 
tain pure emission line fluxes. Contin uum subtracted Hg an d Hj6 
images of IZw 18 are e.g. shown in ICannon et al. I (I2002h . We 
measured a total Ha flux (uncorrected for reddening) of 3.28 x 
10"'^ erg S-' cm-^ within a circula r aperture of 1 .5 "ra dius, 
in agreement with v alues found by Ide Mello et al.l (Il998h and 
ICannon et al](l2002l) . 

Extinction The extinction map of IZw 18 is created using 
the ratio between Ha and H/? images. In the absence of dust 
extinction the theoretical value of the Balmer ratio is known 
to be loosely sensiti ve to temperature and density. Following 
ICannon et all (|2002|) we adopted a value of Viainp^ 2.76. 
Potential sources of error on the expected value, such collisional 
excitation of Hi or underlying stellar absorption, are also ad- 
dressed in this paper An E(B-V) map was thus constructed using 
the following relation: 



E(B - V) Ha I Hp - 



2.5 X log{2.76/Rohs) 
k(Aa) - k(Ap) 



(1) 



where Roi,s = fna/fn^ is the absolute observed flux ratio, and 
k(Aa), k(Aj3) are the extinction curve values at Ha and H/? wave- 
leng ths respec tively. We adopt k(Aa) - kiAp) - -1.08, from 
.Cardelh et al.l (,1989 ). which is appropriate for resolved galax- 
ies and diffuse interstellar regions. Finally, the e xtinction value 
was co rrected for galactic contribution following ISchlegel et al.l 
(Il998h . which accounts for 0.032 mag. The resulting extinction 
map is shown in Fig.H] 

Overall the extinc tion in IZw 18 an d in its NW region is 
known to be very low. lMas Hessd (Il990h found that the Balmer 
decrement of the whol e NW re g ion is consistent with no ex- 
tinction. However, Duf our et alj (Il988h reported an extinction 
of E{B - V) ~ 0.17 in their 2.5"x 6" slit. Ground-based 
spectroscopic observations re vealed typical values ranging from 
E(B-V) ~ 0.03 up to 0.2 (e.g.lVflchez & Iglesias-Paramol 19981: 
llzotov et alJ[T997l:lMartinlll996h . The main reasons of such dis- 
crepancies may be differences in the aperture size and the loca- 
tion of the slits on the galaxy, as the dust does not seem to be 
homogeneously distributed in IZw 18 (see Fig.|4]i. In the present 
work, we define a circular aperture (3.2" radius) centered on the 
NW region. We exclude ( 1 " circular aperture) the central region, 
where Ha and Wfi emission are much weaker and Balmer ra- 
tio gives unreasonably low values, from our measurement. The 
mean color excess derived in this way is E (B -V) ^ 0.042. This 
value agrees with the determinations by ICannon et al.l (|2002|) 
{E{B -V) = - 0.09) obtai ned in different parts of the NW 
region, and with Peguignoj ( l2008l) . We also find no extinc- 
tion when the cen tral region is not excluded, in agreement with 
iMas Hess3 (11990'). Subsequently we will adopt an average value 
of E{B - y) = 0.05 for the NW region. 

Intrinsic Lya emission From the Ha image, and using the 
extinction map, we have produced a theoretical Lya emission 
map (Fig. |4|. It is determined using f{Lya) = 8.7 x f{Ha) x 
]^q(i. 048x£(b-v) )^ where we assume a case B recombination the- 
ory i Brockleh urst 1971) and the extinction law previously men- 
tioned. Naturally, given the small extinction corrections, the re- 
sulting map of intrinsic Lya emission shows a very similar mor- 
phology as the Ha map. Schematically, the NW region is sur- 
rounded by a Lya emission shell including in particular one 
bright knot (NWl in the notation of Cannon et al. 2002)). 

Finally, the spatial variation of the UV continuum, and the 
intrinsic Lya equivalent width are shown in Fig.|5] The UV im- 
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Fig. 4. IZw 18 imaging. Left: Intrinsic Lya emission map in a logarithmic scale. It has been obtained by correcting the observed Ha 
image with the extinction map and assuming case B recombination theory (see text for details). The result has been smoothed using 
a median filter (width=5). Right: shows the extinction map, in linear scale, obtained from the Balmer decrement Ha/H/3 then median 
filtered (width=5). Inverted color scale is used showing higher emission and dust content in darker color. The NW integration box 
(cf. Fig.[Tli is also shown in dashed line. The size of the field of view is about 13"x 15" and the orientation is the same as in Fig.[T] 



age has been coiTected for reddening, scaled to Lya wavelength 
using a UV slope of/? 2, and used together with the intrin- 
sic Lya map to construct the theoretical EW^ya map. Over the 
entire NW region we obtain EW^ya ~ 50 A, compatible with ex- 
pectations for a young starburst. However, as shown in Fig.|5] we 
observe very high values around the UV-bright central region. 

Taken together the observations of strong Lya absorption 
across the entire NW region despite the presence of intrinsic 
strong Lya emission and a very low amount of extinction clearly 
call for a physically consistent explanation of these apparent 
contradicting phenomena. 

2.3. Other observational constraints 

A mean velocity offset, Ay(em - abs), between the systemic ve- 
locity, measured from the optical line s, and metalUc abs orption 
lines of O i and Si ii, was measured bv lKunth et alj(ll998l) in the 
small GHRS aperture centered on the NW region. They found 
Av(em - abs) ~ 25 km s"', indicating that the neutral gas is 
mostly static with respect to the central Hii region. Recent FUSE 
observations including other ISM absorption lines confirm the 
absence of an outflow in IZw 18 on a large apert ure including 
by far all the UV emitting regions of this galaxy; iGrimes et al.l 
(|2008|) measure velocity shifts between ~ and 40 km s~' with 
a mean offset of 8 km s"' . 

The Doppler parameter b describes the thermal motion of hy- 
drogen atoms. Hi velocity dispersion observed by Ivan Zee et al.l 
(119981) is about 12 - 14 km s which translates to ^ 17 - 20 
k m s~'. A slightly hi gher value {b - 27 km s"') was quoted 



by iKunth et al.! (1 19941) from their VLA observations. Given the 
very damped profile of the Lyo- absorption, variations within this 
range of values does not affect the model fit. 

The Full Width at Half Maximum (FWHM) of t he Lyg emis- 
sion line can be constrained using FWHMCHo-). jPufour et al.l 
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Fig. 5. ID emission profile of NW region. The different profiles 
are obtained by adding the flux along the slit (spatial direction) 
in a same aperture used for the extraction of the spectra (cf. Fig. 
|2]i centered on the NW region. Then the ID profile corresponds 
to the direction along the dispersion axis. One pixel corresponds 
to ~ 0.1". The dark solid curve is the 1500 A flux represented 
on the left y-axis, and blue dashed line the intrinsic Lya equiva- 
lent width represented in logarithmic scale on the right axis. The 
position of the seven adjacent slits are also marked, showing the 
spatial variation of the emission between the different slits. 



(119881) found FWHM(Ha) ~ 6. 1 A (280 km s ' ) from their spec- 
trophotometry observations but with a resolution of 275 km s"' . 
Observations with a better resolutio n (R ~ 1 1 km s ' FWHM) 
indicates FWHM(Ha) ~ 150 km s"' (iMartinll 19961) . This is con- 
sistent with a relatively naiTow emission line and we will adopt 
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FWHM(Lya) =100 km s although our results are insensitive 
to the differences found in the observations. 



3. Explaining the Lya absorption in IZw 18 

3.1. General considerations 






L 





Fig. 6. Sketch showing geometrical effects on the Lya profile 
shape for a point source behind a homogeneous dust-free slab 
emitting pure UV continuum radiation (flat continuum around 
Lya). The output spectra represent the observed profiles in dif- 
ferent regions and in the observer's line of sight perpendicular to 
the slab. The reflected spectrum by the slab is also represented. 



To transform the intrinsic Lya emission (emitted in the Hii 
region) to a pure absorption profile can in principle only be 
achieved in two ways: 7 j by true destruction of Lya photons (by 
dust or possibly by conversion to two-photon continuous emis- 
sion in the ionised region), or 2 j by geometrical effects leading 
to the scattering of Lya photons out of the line of sight, or by a 
combination of both. 

E xamples of fine profiles due to dust abso rption are shown 
e.g. in[\^rhamme et alj ( l2006i|200 8^) and S chaerer & Verhammd 
( |2008|) . Effect 2) is illustrated in Fig. |6l showing how for ex- 
ample even a dust-free slab produces an absorption (Voigt) pro- 
file along the central line of sight from a point-like background 
source. This geometrical situation also corresponds to the "clas- 
sical" case of damped Lya systems (DLA) in front of distant 
quasars, or other Lya forest observations. If the scattering fore- 
ground layer was truly dust-free, it is clear that the photons are 
conserved; hence the photons scattered away from line center 
(causing the apparent absorption line) must emerge somewhere. 
In a static configuration, radiation transfer effects redistribute the 
photon s into the wings , leading to a symmetric double peak Lya 
profile (iNeufeldl 19901) . as sketched in Fig.|6]for the distant, non- 
central lines of sight. Adding dust to effect 2), i.e. combining 1) 
and 2), will reduce the strength of the scattered component and 
further increase the depth of the central absorption profile. 

Using radiation transfer models we will now examine 
whether these effects can quantitatively explain the observations 
of IZw 18, and which of these effects is dominant. 



3.2. Lya and UV continuum radiation transfer modeiing 

3.2.1 . MCLya code and input parameters 

We use an improved vers ion of the Mon te Carl o radiation 
transfer code MCLya of IVerhamme et al.l ([2006) including 
the detailed physics of Lya line and UV continuum trans- 
fer, dust scattering, and dust absorption for arbitrary 3D ge- 
ometries and velocit y fields. The following improvements have 
been included (see iHayes et al] 2009, for more details): an- 
gular redistribution functions taking qua ntum mechanical re- 
sults for Lya into account (cf. .Diikstra & Loebll2008t IStenflol 
Il980h . freq uency changes of Lya photons due to the recoil 
effect (e.g. IZheng & Miralda-Escudd l2002h . the presence of 
deute rium (assuming a c anonical abundance of D/H = 3 x 
10"^ iDiikstra et alj 120061) . and anisotropic dust scattering us- 
ing the He nyey-Greenstein phase function (using parameters 
adopted in IWitt & GordonI (|2000|) ). Furthermore a relatively 
minor bug in the angular redistribution of Lya photons has 
been fixed, and the code has been parallelized for efficient use 
on supercomputers. For the physical conditions in the simu- 
lations used for the present paper, these improvements lead 
on ly to minor changes with respe ct to the MCLya version use d 
by ISchaerer & Verhammd ([2008") and IVerhamme et al.' (2008^ 
More details on the code upgrade will be given in IHayes et all 
(l2009l) . 

For simplicity, and given the available observational con- 
straints, all simulations carried out subsequently assume a ho- 
mogeneous and co-spatial distribution of neutral hydrogen and 
dust with a constant density and temperature. The corresponding 
microscopic Hi velocity distribution is described by the Doppler 
parameter b. The remaining input parameters of the code are the 
Hi geometry and velocity field, the spatial location and distribu- 
tion of the UV continuum and line emission source(s), and the 
dust-to-gas ratio. 

We consider the following Hi geometries: spherically sym- 
metric shells with a central source, and plane parallel slabs with 
a background or internal source (including different source ge- 
ometries). These cases are described by 3 additional parame- 
ters: (i) the expansion velocity of the shell, Vexp, or the veloc- 
ity of the slab with respect to the source, (ii) the Hi column 
density towards the source, A^ei, and (Hi) the dust absorption 
optical de pth Tg which expresses the dust-to-gas ratio. As dis- 
cussed bv lVerhamme et alT (l2006l) Ta is related to the usual color 
excess E{B - VO by E(B - V) ~ (0.06...0.1 1)t„; we assume 
E(B - V) = 0.1 Ta for convenience. In short, for a given ge- 
ometry we have 4 parameters (b, Vexp, A^hi, b - 2Q km s"' 
and v 'exp ~ km s"' are constrained by the observations (see 
Sect. 12.31 ), Ta is varied between (no dust) and 0.5, the maxi- 
mum allowed by the observations (Sect. 12.31 ). and A^hi is varied 
to reproduce the observed Lya line profile. 

For each parameter set a full Monte Carlo simulation is run 
allowing for sufficient statistics to compute both integrated and 
spatially resolved spectra in the Lya region. The radiation trans- 
fer calculations cover a sufficiently broad spectral range (here 
typically from - 1 0000 to -i- 1 0000 km s" ' ) necessary to reach the 
continuu m for the highest column density simulations. As de- 
scribed in lVerhamme et alj (12006) our MC simulations are com- 
puted for a flat input spectrum, keeping track of the necessary in- 
formation to recompute a posteriori simulations for arbitrary in- 
put spectra. For the Lya fits we assume an input spectrum given 
by a flat (stellar) continuum plus the Lya line, described by a 
Gaussian with variable equivalent width EW^/a and full width at 
half maximum FWHM{]jya). EWi^ya is kept free, although con- 
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straints are available from our theoretical (intrinsic) Lya map; 
a FWHM - 100 km s ' is assumed as for Hot, although our 
results are basically independent of its exact value. Other con- 
tinua, such as syntheti c high resolution starburst spectra from 
ISchaerer & Verhammd (l2008h . can also be used. 

3.2.2. Shell models 

To consider a simple geometry to understand the observed Lya 
absorption of IZw 18, we examine predictions for the integrated 
spectrum of a spherical shell with/without dust. In this case no 
"loss" of photons by spatial diffusion is allowed. Hence to trans- 
form intrinsic Lyor emission into an absorption profile requires 
absorption by dust. We will now examine whether spatially in- 
tegrated shells can recover the observed profile for reasonable 
amounts of dust and reasonable Hi columns. 



— 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 J. 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 

■ observed spectrum - 
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Fig. 7. Comparison of the observed and fitted Lya profile of IZw 
18 assuming a spherical shell model. The observed STIS spec- 
trum (from MHOS) is represented by the dark line. The model fit, 
using a shell geometry, is plotted in red dashed line and the pa- 
rameters used are: A^hi = 6.5 x 10^'cm"^, r = 0.5, v=0 km s"', 
b = 20 km s The blue dotted line represents the input spec- 
trum of the simulation. It consists of a flat UV continuum plus 
a Gaussian Lya emission line with FWHM-\QQ km s The 
intrinsic equivalent width adopted is EW(Lya)-()Q A. 



net absorption is only caused by dust absorption, since, in the 
present case, we observe all the scattered photons escaping from 
the shell, without any hne-of-sight effect. 

Influence of V^^p: We adopted a static shell in our model to 
fit the Lya absorption profile. As discussed in IVerhamme et al.l 
(l2006h . for increasing Vexp more Lya photons will escape from 
the red part of the line, because Lya photons are seen already 
redshifted by hydrogen atoms. However, since the high column 
density reduces the escape probability, we can vary the expan- 
sion velocity in a certain range without affecting the quality of 
the fit. The highest velocity allowed is around 50 km s"', which 
already exceeds the observed outflow velocity of IZw 18. 

Influence of ^WLy^: We can show that EW \^a close to the m ax- 
imum value expected by synthesis models ( ISchaererll2003h . for 
normal IMF populations, are allowed in our model to fit IZw 18 
profile. Indeed, when we use an input Lya line with EWij^a - 
200 A the fit remains globally unchanged. Again, in spite the 
large damped absorption, high intrinsic EW]_ya is not excluded 
by radiation transfer simulations because of the high Hi column 
density. 

Other solutions: We need to invoke a relatively high column 
density (A^hi - 6.5 10^' cm"^) to obtain a good fit of the ab- 
sorption wings. On the other hand, all photons which, in reality, 
will scatter away from the observer's line of sight, are recov- 
ered in our simulation, since we integrate over the entire surface 
of the shell. Relaxing this assumption, i.e. considering different 
geometries, would in particular also allow us to lower A^hi- 

The solution proposed here to fit IZw 1 8 profile is not unique, 
and different combinations of A^hi and t„ can reproduce the ab- 
sorption. For instance, the use of a higher value for the extinction 
(ja- 1) and a lower Hi column density (A^hi= 5x 10^' cm"^) pro- 
duce the same fit quality. Overall, this somewhat academic case 
of a shell model for IZw 18 serves to show that even low dust 
quantities may suffice to transform Lya emission into a broad 
absorption profile, provided a sufficiently high column density 
and/or sufficiently low outflow velocity, as also discussed in 
Sect.m In any case, the radio observations of IZw 18 show very 
clearly a large spatial extension of Hi compared to the size of 
the NW region (and to that of the spectroscopic apertures). The 
effect of such geometries on Lya will be addressed now. 



Adopting an average extinction of E{B - V) ~ 0.05 (i.e. 
Tfl = 0.5) and b - 20 km s"' we have computed several static 
shell models with varying A^hi- As shown in Fig.|7] the predicted 
profile agrees well with the observations for A^hi = 6.5 x 10^' 
cm"^ and for an input spectrum with a Lya line equivalent width 
^^(Lya) = 60 A. Note that we do not fit the absorptions in 
the blue wing, attr ibuted to Si n AA 1193.3,119 4,5 and Si iii 
AA 1206.5,1207.5 (ISchaerer & V erhamme' '2008*) and possibly 
to Galactic and intergalactic Hi absorption, since these are not 
taken into account in our model. 

The reason for the resulting broad damped Lya absorption 
is as follows: Due to the high Hi column density, even a small 
amount of dust destroys almost all photons in and around the 
Lya line center. Scattering on hydrogen atoms with such a high 
column density, greatly increases the mean path of Lya photons, 
and hence the probability to be absorbed by dust. Therefore, the 



3.2.3. Extended geometries and line of sight effects 

The galaxy spectrum we observe in reality, could deviate signifi- 
cantly from the simple homogeneous shell model presented here, 
since the source is spatially resolved and the spectrum is not in- 
tegrated over the whole shell surface. Furthermore the spectrum 
can depend on viewing angle and on the geometry of the ISM. 

We show in Fig. [8] that IZw 18 absorption can be well ad- 
justed with lower Hi column density than required for the shell 
model and without any dust (E(B - V) - 0). This result is 
achieved by taking a slab geometry with a static gas and ap- 
plying sight-line selection criterium, where only photons in the 
observer's direction are collected. Then the absorption is caused, 
not by dust destruction, but by diffusion of the photons out of the 
observer's direction. Strictly speaking, no photon is destroyed. 
This demonstrates even better that Lya absorption can be ob- 
served in dust-free galaxies (cf. Fig.|6]for a schematic overview). 
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Fig. 8. Lya absorption fitting II. The observed spectrum (STIS 
MH03) is represented by the dark line. A slab geometry is used 
for the model spectrum and only photons in the observer's sight- 
line are collected. It is plotted in green dashed line and the pa- 
rameters used are: A^hi = 3 x 10^'cm"^, t = (no dust), v=0 
km s b = 20 km s~'. The blue dotted line represents the input 
spectrum of the simulation. It consists of a flat UV continuum 
plus a Gaussian Lya emission line with FWHM-IOQ km s"'. 
The intrinsic equivalent width adopted is £'W(LyQ;)=60 A. 



Only a nearly static neutral ISM is required, with A^hi ~ 3 x 
10^'cm"^ in this case. 

In this case, we expect to recover the diffused photons in 
other directions and/or further from the source. On the other 
hand, in presence of dust, this diffuse part would be attenuated or 
suppressed. For example, for models with homogen eous gas and 
dust d istributions, our Lya transfer simulations (see lHaves et al.l 
|2009|) predict already quite low escape fractions for Lya line 
photons, with /esc of the order of typically 5-10 % for column 
densities Am ^ 10^' cm"^, dust optical depths Tq = 0.2, and 
low expansion velocities (vexp ^ 50 km s"'). Much lower escape 
fractions (/esc ~ 10"^ are predicted for larger amounts of 
dust, such as for the average value adopted for the NW region. 
Therefore we expect relatively small amounts of diffuse emis- 
sion from Lya line photons. 

3.2.4. Spatial variations of Lya profile 

In Sect. 12. II (Fig. |2Ti we have shown that the Lya profile shows 
spatial variations between the different STIS slits. We now 
demonstrate, that given the observational constraints, the Lya 
radiation transport explain fairly well these variations. We con- 
sider a large, static and uniform cloud of Hi (Ahi = 3 x 10^' 
cm"^, b - 2Q km s"', Ta- 0.5) represented by a slab geom- 
etry, covering the NW star-forming region. We then simulate 
the observed spatial variations of the emission strength by using 
weighted point sources located in front of the Hi slab, emitting 
a flat UV continuum, as input to our radiation transfer code fol- 
lowing the observed UV profile of Fig. |5] The addition of Lya 
line emission will be discussed below. 

The result of this simulation is shown in Fig. |9] The output 
spectrum is what an observer would see when his line-of-sight 
is perpendicular to the slab surface. At the center (1 150 < pixel 
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Fig. 9. Predicted spatial variations of the Lya profile. The simu- 
lation consists of an homogeneous Hi slab with Ahi = 3 x 10^' 
cm"^ and = 0.5 illuminated with an isotropic, extended, UV 
source emitting a flat continuum. The emission strength is spa- 
tially varying from the center to the edge to reproduce the ob- 
served UV surface brightness of IZw 18. An observer line-of- 
sight perpendicular to the cube is chosen {9-0). The different 
plotted profiles correspond then to different regions at the sur- 
face of the slab and are marked in units of pixels corresponding 
to the position of the slits in Fig. |5] The observed profiles (cf. 
Fig.|2]l are plotted in dotted lines. The blue wing of the profiles 
is not well fitted because it is affected by geocoronal emission 
and Hi Galactic absorption which the model does not account 
for. 





V 



Fig. 10. Geometrical effects on the Lya profile shape: 
Homogeneous, static, and dust-free slab of neutral gas illumi- 
nated by a series of isotropic point sources emitting UV contin- 
uum radiation centered on Lya wavelength. An extended source 
is simulated in this case with a varying emission strength, sym- 
bolized by the size of the individual point sources. The output 
spectra represent the observed profiles in different regions and 
in the observer's line of sight perpendicular to the slab. 



< 1 155), in the direction of the brightest source, we observe the 
strongest (largest) Lya absorption profile. The profile proves in- 
creasingly narrower as one moves away from the center, what 
reproduces the trend observed in IZw 18. The double -peak con- 
tribution, characteristic of diffused photons, can even be seen in 
the peripheral region (1 135 < pixel < 1 140). 
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To understand these results, let us decipher the different con- 
tributions in the simulation. Figure |6] depicts the situation for 
this purpose. It shows the observed spectra in a simulation using 
a point source and isotropic emission behind a uniform slab of 
neutral gas. Observing the slab face-on, toward the source, we 
obtain an absorption profile. Only photons far from the line cen- 
ter are transmitted directly, forming the "continuum". Photons in 
the line center are resonantly absorbed and reemitted, diffusing 
in frequency and in space, and leading to the lack of emission at 
and around the line core. These photons will be collected if we 
look at the cube at a position far from the source. A double-peak 
profile is then observed consisting of the diffused photons and 
the absence of photons that would have escape directly, without 
scattering, in this direction. 

Figure [To] now shows a combination of these single sources 
but with different intensities, illustrating the extended source 
simulation of Fig. |9] As for the single source, the spectrum of 
the central region shows a typical damped absorption. At the po- 
sitions of the fainter peripheral sources two contributions lead to 
a narrower absorption profile: a) The transmitted flux is fainter 
than in the central region, and b) photons that have diffused from 
the brighter sources to escape further (double-peak emission), 
contributing to "fill the wings". In the central region the direct 
transmission is stronger and the diffuse part is weaker. In this 
way spatial variations of the UV continuum combined with the 
resonant transport effects of Lya radiation, can explain qualita- 
tively the observed Lya profile variations in IZw 18. The ob- 
served profile in slit 7 (Fig. |2]i, may even show a hint of the 
predicted double peak profile in its red wing, although the S/N 
is quite low in this region. 

Now we discuss the effect of adding Lya line emission on top 
of the UV continuum emission. Naively one could expect a very 
different behavior given the very large Lya equivalent width of 
the source in the peripheral parts of the NW region (cf. Fig. |5]l. 
However, the final spectrum remains unchanged despite the high 
Lya equivalent width used. It appears that the photons emitted 
at the core of the line are either destroyed by dust (Ta= 0.5 here) 
or backscattered, and only photons with higher frequency shift 
diffuse and contribute to the double peak emission. Therefore, 
increasing EWi,ya has no incidence on the output spectrum since 
with FWHM(LyQ')=100 km s all photons are emitted close to 
the center This is easily confirmed by looking at the reflected 
spectrum (cf. Fig. |6]l which increases with higher EW^ya- We 
need to use unreasonably high FWHM(LyQ') (^ 1000 km s"') to 
affect our result and see the double peak contribution increas- 
ing (in the profile wings). This implies in particular also that our 
model predictions are insensitive to the observed spatial varia- 
tions of EWi^ya (cf. Fig. Is). 

In short, we conclude that the observed variations of the Lya 
profile across the NW region can be understood by a combina- 
tion of the line of sight effects discussed earlier and by radiation 
transfer effects related to an extended source. 

3.3. Discussion 

For the first attempt to reproduce the damped absorption profile 
of IZw 18, we used a simple expanding shell model (Sect. [TZ2l l. 
If line-of-sight arguments could not be invoked, we would need 
a relatively high column density (A^hi = 6.5 x 10^' cm~^) and a 
minimum amount of dust {E{B -V) ~ 0.05), which in this case is 
the only way to lose Lya photons. However, when we spatially 
select photons in the observer's sight-line, we showed (Fig. [8]l 
that one may observe Lya in absorption even without any dust 



(E(B -V) - 0). These conclusions also hold for the SB region of 
IZw 18 for which the integrated spectrum show a slightly larger 
Ly a absorption (a Voigt fit yields A^hi ~ 4 x 10^' cm"^). 

iMartinI d 19961) found evidences of supergiant shell in IZw 
18 expanding at a speed of 35-60 km s The geometry pro- 
posed is a bipolar shell seen almost perpendicularly to its main 
axis (cf. their Fig. 4). This configuration is comparable to the 
shell geometry adopted here (Sect. [JX2] i given the negligible ef- 
fects of such small expansion velocities on our model spectrum. 
However, the output spectrum of the shell model would be sig- 
nificantly affected if the Hi coverage is inhomogeneous and low 
column densities are observed in some sight-lines, which is still 
unclear here. For same reasons (low velocity and large Hi cov- 
erage), applying our extended geometry scenario (Sect. |3.2.3t 
to this configuration would yield same results, as our sight-line 
selection is still compatible with this ISM morphology. 

From our shell model we derived an Hi column density of 
6.5 X 10^'cm"^ which is higher than indep endant measurement s 
of Am- 2.6x10^' cm^^ fromradio databv lvan Zee et al.l(ll998l) . 
Since, radio observations measure the total Hi content, our 
model would imply, for symetrical reasons, a value twice higher 
than given. However, with a typical beam size of 5", the radio 
observations are not able to resolve the potential subparsec-size 
Hi clumps in the NW region, and the smoothing effects could 
lead easily to an under-estimation by a factor of 2 -4. Therefore, 
in absence of higher resolution observations, we can not rule out 
the supershell geometry. 

The second model adopted led to a good fit with Ahi= 3 x 
10^' and without dust comparable to observational constraints, 
and also explained the spatial variations of the absorption profile. 
It is therefore more likley that if the emission region is embed- 
ded in an Hi region, the geometry would be not symetric, with 
a higher column density in the front and/or ionised holes in the 
backside. Finally, it is worth noting that the geometry proposed 
by Martin (1996) is not an embedded-like source and the ex- 
panding shell is bipolar and asymmetric with an a xis inclined by 
/ ~ 10° to our line of sight (see their Fig. 4), while Ivan Zee et all 
(119981) find a higher incUnation of / ~ 55°. 

4. Comparison of IZw 18 with other nearby and 
high-z starbursts 

We have just shown how with a low extinction or even no dust at 
all it is possible to explain by radiation transfer and geometrical 
effects the transformation of a strong intrinsic Lya emission into 
the broad Lya absorption profile observed in IZw 18. We need 
now to understand whether this galaxy is unique or representa- 
tive of a certain class of objects and what our results imply for 
other studies, including in particular Lya observations of high-z 
objects. 

4.1. Comparison witli local starbursts 

Four of the eight Hii galaxies observed with GHRS/HST by 
iKunth et"an [l998) show broad Lya absorption profiles: II Zw 
70, Mrk 36, SBS 0335-052, and IZw 18 studied here. As already 
noted by these authors, these objects clearly differ from those 
with Lya in emission by very low velocity shifts between the 
interstellar absorption lines and the systemic velocit}Q. An es- 

- Two of the objects with Lya in absorption, IZw 18 and SBS 0335- 
052, have also b een observed with FUSE, with a large aperture. The 
measurements of lGrimes et al.l 12008) confirm the earlier finding of low 
velocity shifts, now also on a much larger aperture. 
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sentially static ISM appears therefore as one of the main factors 
leading to Lyor absorption, as already concluded by these authors 
and as supported by our radiation transfer modeling. 

Furthermore, among the Lya absorbers, SBS 0335-052 and 
II Zw 70 show clearly higher extinction, with E(B - V) = 0.18 
and 0. 15 respectively (less than 0.02 for Mrk 36. llzotov & Thuaiil 
(Il99 8l)). Hence the ISM properties of these objects should ful- 
fill the same conditions, which have allowed us to explain the 
Lya absorption of IZw 18, and dust destruction of Lya photons 
should be equally or more important. Although very likely, we 
cannot fully prove this until Hi column density measurements 
from the radio are available for all of them. For SBS 03 35-052 
A^Hi re aches high values, up to ^ 9.4 x 10^" cm"^ ( Pus tilnik et al.l 
1200 lb . Similarly, Mrk 36 shows a high colu mn density peak up to 
2.4 X 10^' cm"^ (iBravo-Alfaro et al.ll2004:) . Of course, depend- 
ing on the efliciency of dust destruction, some spatial regions 
with diffiise Lya emission may be expected; however, this is not 
necessarily the case. For example, for SBS 0335-052 we know 
that Lya absorption is observed over a large area, show ing that 
absorption by dust must be important ('Atek et al."20 08h. 

The other half of the HST sample of Kunth etlD ( Il998l) 
shows Lya profiles in emission and varying amounts of dust 
(E(B - V) ranging from ~ 0.02 to 0.18). As akeady mentioned by 
these authors, the main difference with the other part of the sam- 
ple showing Lya absorption appears to be the clear signature of 
ISM outflows in the former A continuity of ISM velocities be- 
twe en "static" and "outflowing" is ex pected and observed (see 
e.g. iMartinI 120051: iGrimes et al] 12008 ). mosdy correlated with 
galaxy luminosity, stellar mass, and SFR. A more detailed anal- 
ysis of the full sample o f nearby starbursts observed in Lya will 
be presented elsewhere (lAtek et al...2009.) . 

4.2. Comparison with distant galaxies 

Compared to distant galaxies it is clear that IZw 18 and SBS 
0335-052, or at least the regions of these objects showing intense 
star formation, show very high Hi column density. For example, 
with A^Hi ~ (0.9 - 3) X 10^' cm"^, these two regions would corre- 
spond to the hig h A^hi tail of all DLA sy stems found in the SDSS 
DR3 survey (cf. iProchaska etaLll2005b 

Also, few high-z starbursts with Lya absorption as broad 
as in IZw 18 an d SBS 0335-052 are known. While ~ 25% 
of the LBGs of Sha plev et al.l (1200 3*) show Lya absorption, 
their stacked spectrum shows a narrower absorption profile. 
Among the broadest Lya profiles of z ^ 3 LBGs are the two 
lensed galaxies MS 1512-cB58 and FORJ0332-3557, whose 
absorption profiles corresponds to ~ ( 0.7 - 2.5) X 10^' cm"^ 
dPettini et al.ll2000t ICabanac et af 2008). 

However, LBGs in general and these two objects in partic- 
ular diff'er in many properties compared to IZw 18. The ob- 
jects with strong Lya absorption show significant dust extinction 
(E^iB -V) ~ 0.169 + 0.006, where E^{B - V) is the color ex- 
cess determined from stellar light), and high Star Formation Rate 
(SFR) (dust-corrected SFR ~ 52 + 5 M© yr_i). Furthermore out- 
flows with significant velocities (vgxp ~ 100 - 300 km s"') are 
generally observed in LBGs. In comparison, IZw 18 is a very 
low luminosity, low SFR object (with a UV luminosity lower 
than that of LBGs by 2-3 or ders of magnitude, SFR(UV) ~ 0.3 
Mq yr ' dGrimes et al.ll2008l) with a low extinction (E(B - V) < 
0.05), which s hows a static ISM. 

For LBGs ISchaerer & Verhammd (l2008h have shown with 
radiation transfer models that the absorption profile of MS 1512- 
cB58 is due to the relatively large amount of dust and the high 
column density; with the observed ISM conditions this suffices 



to transform intrinsic Lya emission expected from the ongoing 
starburst to broad Lya absorption, despite the relatively large 
outflow velocity (vexp ~ 220 km s"'). In IZw 18 geometrical 
effects or a static high A^hi ISM with small amounts of dust are 
sufficient to do a similar "transformation". 

In short, we suggest schematically the following two expla- 
nations for Lya absorption in nearby and distant starbursts: 1) 
On average the cold ISM (relevant for Lya transfer) of LBGs 
shows the geometry of a spherically expanding shell with rel- 
atively large velocities and smal l variations (Vexp ~ 10 0-300 
km s"^) (cf. ISh apley et al. 2003; ISchaerer & Verhammd 120081; 
IVerhamme et al.l 120081) . In such cases the main factors deter- 
mining the escape fraction of Lya ph otons are A^hi and Tg , 
as shown by radi ation transfer models (IVerhamme et al J [20081 ; 
iHaves et al.ll2009l) . and significant amounts of dust are required 
to obtain broad Lya absorption profiles. 2) In nearby galaxies, 
small amounts of dust in a static/low velocity ISM with a high 
Hi column density suffice to create Lya absorption. Furthermore, 
the occurrence of Lya absorption is most probably metallicity 
independent, at least to first order In addition, geometrical ef- 
fects due to small apertures may also increase the observed Lya 
absorption. 

The distinction between groups 1) and 2) is most likely sim- 
ply due to the outflow properties, i.e. the wind velocity, which is 
known to increase with SFR, galaxy mass, and speci fic star for- 
matio n rate (e.g. Martin 2 005; Rupke et al, 2005; S chwartz et al.l 
20061; l&imes et all2008l) . Qualitatively this increase of the out- 
flow velocity behavior is understood by increasing mechanical 
feedback on the ISM related to stronger SF activity (SFR) in 
galaxies with increasing mass or luminosity. At the low lumi- 
nosity (SFR) end, feedback appears to be insufficient to "ignite" 
outflows, hence the nearly static ISM in IZw 18 and alike ob- 
jects. What ultimately settles the ISM geometry, A^hi and dust to 
gas ratio, and hence assures in particular a high Hi column den- 
sity in IZw 18 and other local objects remains to be explained. 

Clearly, the observed trends and diversity need to be exam- 
ined further both qualitatively and quantitatively. This will be the 
scope of subsequent publications. 

5. Summary and conclusion 

Archival HST/STIS UV spectroscopy and imaging, and 
HST/WFPC2 optical imaging data of the nearby star forming 
galaxy IZw 18 were obtained. We have appli ed the 3D Monte 
Carlo Lya radiative transfer code MCLya (IVerhamme et al.l 
12006 ^ to explain quantitatively the intriguing Lya absorption in 
this galaxy and the apparent spatial variation of the Lya profile. 
Then, using the example of IZw 18, we have discussed under 
which physical conditions one observes Lya in emission or ab- 
sorption both in nearby or high-z galaxies. Our main results can 
be summarized as follows: 

• We first examined the predictions of a spherical shell model 
to reproduce the integrated spectrum of the NW region of 
IZw 18. This model described a static shell of Hi mixed with 
dust grains, surrounding a central point source emitting UV 
continuum plus a Lya emission line. Adopting dust extinc- 
tion derived from observations {E{B - V) = 0.05) and b -2Q 
km s ', we were able to fit the Lya profile with Hi column 
density of A^hi= 6.5x10^' cm"-. Even with a strong intrinsic 
Lya emission EWi,ya (up to 200 A) a small amount of dust 
is sufficient to cause strong damped Lya emission, since the 
probability to be absorbed by dust is greatly increased by 
the high column density and by a nearly static ISM. In this 
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model, the loss of Lya photons is only possible by means of 
true dust absorption, since we spatially recover all photons. 

• Give n the large spatial e xtension of Hi covering the NW re- 
gion jvan Zee et al.ll 1998) . we also explored the slab geome- 
try of neutral gas in front of the UV source. We have shown 
that considering only emission along the observer's line of 
sight, we can reproduce the strong Lya absorption without 
any dust {E{B - V) - 0). This is achieved by the diffusion 
of the photons out of the observer's sightline provided a suf- 
ficiently high column density (at least A^hi= 3 x 10^' cm"^) 
and a nearly static ISM configuration. 

• We have observed spatial variations in the Lyff profile shape 
in the different STIS slits. From observations we have con- 
structed the ID profile variations of the UV continuum and 
EWi,ya across the NW region. Despite strong UV emission 
at the center and high EWi,ya around, the Lya profile is still 
in absorption in all the NW region and proves narrower to- 
ward the peripheral region. To understand these variations 
we have simulated an extended source with a spatially vary- 
ing UV emission strength by following the ID spatial pro- 
file, in front of a slab of neutral gas with A^hi= 3 x 10^' cm"- 
and E(B - V) - 0.05. Then, by observing the slab at differ- 
ent distances from the center, we have been able to repro- 
duce the observed spatial variations. We have demonstrated 
qualitatively that this is due to radiative transfer effects, in 
particular, to the diffusion of Lya photons and to the spa- 
tial variation of the UV continuum source. Finally, adding 
a Lya recombination line to the source simulating the ob- 
served EWi^ya profile, does not affect the final spectrum. 

Other nearby galaxies with intense star formation, such as II 
Zw 70, Mrk 36 and SBS 335-052, also show strong Lya ab- 
sorption (iKunth et al.lll998l) . They appear to show a nearly static 
cold ISM and more or less significant amounts of dust. At least 
two of these objects exhibit very high Hi column density (up 
to 2.4 X 10^' cm"^). We suggest that Lya absorption in these 
objects is due to the same reasons as for IZw 18: a very large 
number of scatterings in static, high column density gas leading 
to an efficient destruction of Lya photons by even small amounts 
of dust. Furthermore line-of-sight effects can also be responsible 
for or increase Lya absorption further. 

The distinction between Lya emission and absorption in lo- 
cal starbursts seems to be mainly related to presence or not of 
ISM outflows. Since high-z objects (LBGs, LAEs) show gen- 
erally outflows with high, but relatively similar velocities (with 
bulk velocity typically ~ 100-200 km s '), the transition from 
Lya absorption to emission in these objects is, on the other hand 
mostly determined by the dust content and Hi column density 
(cf. Schaerer & Verhamme 2008; Verhamme et al. 2008). 

These results and the global trends observed between Lya 
strength and profile diversity, and relevant parameters need now 
to be tested with larger samples of gal axies. This is the main 
objective of upcoming publications (eg. I Atek et al.ll2009l) . 
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